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A highly collimated, extremely high velocity outflow in Taurus 
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Abstract. We present the first case of a highly collimated, extremely high velocity bipolar outflow in Taurus. It is powered by 
the low-luminosity (0.4 Lq) source IRAS 04166+2706 and contains gas accelerated up to 50 km s"' with respect to the ambient 
cloud both toward the blue and the red (uncorrected for projection). At the highest velocities, the outflow collimation factor 
exceeds 20, and the gas displays a very high degree of spatial symmetry. This very fast gas presents multiple maxima, and most 
likely arises from the acceleration of ambient material by a time- variable jet-like stellar wind. When scaled for luminosity, the 
outflow parameters of IRAS 04166 are comparable to those of other extremely high velocity outflows like L1448, indicating that 
even the very quiescent star-formation mode of Taurus can produce objects powering very high energy flows (L,„^£./L, > 0.15). 



1. Introduction 

The ejection of supersonic gas in a bipolar outflow is one 
of the first signatures of stellar birth and very likely a re- 
quired element in the physics of star formation. Outflows have 
, been found towards young stellar objects (YSOs) of nearly all 
' masses, and in most cases, their molecular component (traced 
in CO) presents a weak degree of c ollimation and a preva- 
lence of low velocity gas JLadalllQSSL) . A small group of out- 
flows, however, appears dominated by an extremely high ve- 
■ locity (EHV) molecular component having a speed over 20 km 
, s" ' (Mach > 1 00 for gas at 1 K) and a very high d egree of col- 
limation (length/width > 10, see iBachiller Ill996l) . This group 
of highly collimated EHV outflows seems associated with the 
youngest protostars (Class O. Andre et al. -1993.) . and probably 
represents a special ly energetic initial stag e in the evolution of 
a bipolar flow (e.g.. lBontemDS et al. 119961) . 

Although the nearby Tauni s cloud contains a large popula- 
tion of outflows ([Hey^ret al. Ill9 87: Moriartv-Schieven et al. 
Il992t lAndreetal.lll999l) . no Taurus YSO has been found 
so far to power an EHV flow. This is surprising since other 
nearb y clouds like Perseus and Ophiuchus contain EHV 
flows (Bachiller e t anil990t \l99 T: McCa ughrean et al. lfl994l 
lAndre et al. .19901) . and it may at first suggest that the gas con- 
ditions in Taurus prevent the formation of this type of systems. 
In this paper, however, we report the identification of the first 
highly collimated EHV outflow in Taurus, which is associated 
with the low luminosity (0.4 Lq) source IRAS 04166+2706 
(104166 hereafter). The youth of this s ource has been pre- 
viously r ecognized ( Hg^^e^^lJ_|l222l |KgnXQng^^lJ[ll22^ 
Barsonv & Kenv oiil Il992t jChenetaLl Il995[ Lshirlev et al. I 
2000t iMotte & Andre 11200 it IPark & Kenvon 112002' ). but only 
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a poorly collimat ed outflow and weak [S I I] emission have bee n 
reported so far dBontemps et al. Ill996l lOomez et anil997t) . 
In the following sections we present the newly found EHV 
outflow and discuss some implications for outflow and star- 
formation studies. A future paper will describe the shock- 
induced chemistry of this flow (Santiago et al. 2004, in prepa- 
ration). 



2. Observations 

We observed CO(J=2-l) with the 9-pixel HERA receiver array 
on the IRAM 30m telescope in 2004 April. In order to obtain 
flat baselines, we observed in wobbler switching (WSW) mode 
chopping 240" away in azimuth. This procedure contaminated 
the spectra at low velocities due to the extended Taurus emis- 
sion, so we also observed in position switching (PSW) mode 
using a reference more than 1 degree away and free of CO 
emission at the velocities of interest; these observations have 
been used to reconstruct the emission at low velocities. In ad- 
dition, 4 selected positions were observed in CO(l-O) using 
WSW mode to help estimating excitation conditions. All these 
data were taken with the VESPA coiTelator set to velocity res- 
olutions of 0.4 and 0.8 km s' for CO(2-l) and CO(l-O), re- 
spectively. Cross scans of bright continuum sources were used 
to coiTect the telescope pointing, and standard efficiencies were 
used to convert the data into the main beam brightness scale. 
The FWHM of the telescope beam was 22" and 11" for CO(l- 
0) and CO(2-l), respectively. 

We observed NH3(1,1) and NH3(2,2) simultaneously to- 
ward 104166 with the 100m telescope of the MPIR at Effelsberg 
in 1998 October. We used frequency switching mode and the 
AK90 autocorrelator to achieve a velocity resolution of 0.03 
km s '. Cross scans of continuum sources were used to correct 
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Fig. 1. CO(2-l) integrated intensity maps for the EHV, IHV, and SHV regimes (see text). Note the high symmetry between the 
EHV blueshifted gas (blue contours to the NE) and the EHV redshifted gas (red contours to the SW). The EHV and IHV maps 
have been made using WSW spectra, while the SHV map has been made using PSW data. First contour and interval are 1 K 
km s~' in all plots except for the more prominent blue gas in the rightmost (SHV) panel, where they have been increased to 
3 K km s"' to avoid overcrowding. Offsets are referred to the nominal position of IRAS 04166+2706 (q'720()0 = 4*'19™42'!6, 
<5y2000 = 27°13'38"), and Vq - 6.7 km s"' has been derived from NH3(1,1) data. 



the telescope pointing, and an observation of LI 551 was used 
to calibrate the data. The telescope beam FWHM was 40". 

3. The outflow 

3.1. Morphology and kinematics 

Figure 1 presents integrated maps of the CO(2-l) emission to- 
ward 104166 for three velocity ranges. These ranges are offset 
on average from the cloud systemic velocity (= 6.7 km s"') by 
40, 20, and 8 km s"', and will be referred hereafter as the ex- 
tremely high velocity (EHV), intermediate high velocity (IHV), 
and standard high velocity (SHV) regimes, respectively. In all 
maps, the CO emission forms two well-separated lobes that 
have 104166 near their apex, clearly indicating the presence of 
a bipolar outflow powered by this source. At the most extreme 
velocities, the outflow emission is unresolved transversally by 
our 11" beam, so even if we only consider the innermost (con- 
nected) region of emission (220" long), the implied collimation 
factor is at least 20. This factor makes the 104166 system o ne 
of the most coUimated outflows known (see iBachiller 1 1 9961) . 

As Figure 1 shows, each EHV component consists of a con- 
nected lobe 110" (0.07 pc) long near 104166 and additional 
emission at further distance. Both blue and red connected lobes 
break up into peaks almost symmetrically located with respect 
to the IRAS source. Each of these lobes ends in a prominent 
peak, and while the blue lobe contains three additional maxima, 
the red one has two and a hint of a third one in between. These 
presence of symmetric peaks is common in EHV outflows from 
low mass stars, like LI 448 and 103282 (so-called "bullets": 
IBachiller et al. Il99fllll99ll) . and is suggestive of a discontinu- 
ous, even periodic, mass ejection mechanism. Dividing the av- 
erage CO speed of the EHV gas (40 km s ') by the average 
peak separation (20"), we infer an ejection period of 350 yr. 



More than 110" away from the IRAS source, both EHV 
lobes contain additional discrete peaks. The prominent blue 
peak near (85", 145") seems to have a counterpart near 
(-100", -170") (unfortunately not fully mapped due to lack 
of time), and there is a hint of an additional EHV peak in the 
southernmost part of the blue lobe. This peak is probably part 
of a bright IHV component, and may have a weak blue coun- 
terpart at the northernmost part of our map. Overall, the EHV 
emission shows a high degree of blue-red symmetry. 

The IHV emission (middle panel of Fig. 1) shows a more 
diffuse distribution and a lower degree of symmetry (but see en- 
ergetics below). The blue IHV lobe follows the EHV emission 
and deviates slightly to the east at large distance from the IRAS 
source. The red lobe also deviates east at intermediate distances 
but lies again close to the outflow axis by the end of our map, 
where it presents a discrete peak. As mentioned before, there 
is a hint of blue IHV emission at the north end of our map, but 
even if we discount this region (which needs further mapping), 
the IHV extends more than 450", which corresponds to more 
than 0.3 pc at the distance of Taurus (140 pc). Dividing this 
projected distance by the average radial velocity extent of this 
component we derive a dynamical age of 7500 yr, and using 
the highest outflow speed reduces the estimate to 3000 yr. 

Finally, the SHV emission forms two conical lobes with 
the same axis as the EHV gas and relative maxima close 
to the apex. Again there is a slight enhancement toward the 
east, but the bulk of the SHV emission is rather symmet- 
ric and envelopes the EHV emission. This suggests a strati- 
fication of the outflow gas, with the EHV gas moving along 
the axis and the slower SHV material at further distances 
(see lOueth & Guilloteau ' 1999 for a similar case in HH211). 
The overall high collimation of the outflow allows to sepa- 
rate it from the [SII] jet-like feature HH 390 identified by 
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Table 1. Outflow energetics 
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Fig. 2. Top: CO(2-l) spectra toward selected outflow posi- 
tions. Note how the EHV gas appears as a distinct feature in 
the spectra. Bottom: Position-velocity diagram of the CO(2-l) 
emission along the full length of the outflow, also showing the 
distinct nature (and near constant velocity) of the EHV compo- 
nent. First contour and interval are 0.05 K. All data are from 
WSW observations with the emission between V/,,- -3.3 and 
16.7 km s ' reconstructed from PSW observations. Offset ori- 
gin as in Fig. 1. 



iGomez et al. I ( 119971) . which seems powered by a separate 
source almost 2. '5 north of 104166 (J. Bally and J. Walawender, 
private communication). 

3.2. Energetics and central core 

Fig. 2 shows a series of CO(2-l) spectra (top) and a position- 
velocity (PV) diagram along the outflow (bottom) illustrating 
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the three velocity regimes. The slower SHV gas appears as a 

wing in the spectra, while the weaker IHV gas forms a flat, low 

level component more extended toward the south. The EHV 

gas, on the other hand, forms discrete maxima in the spectra, 

like th e EHV features found in o ther highly collimated out- 
j II 1 1 I ° -^ 

flows JBachiller et a l. 1990, 1991). This component seems to 

decrease in velocity slightly along the outflow axis (< 10 km 
s '), and its spectral signature varies significantly with position 
(compare the top two spectra). 

Given the distinct nature of the velocity components (see 
Sect. 4 for further discussion), we estimate their energetics con- 
tent separately. We assume that the CO(2-l) emission is opti- 
cally thin and we do not correct for projection, so our results 
are strict lower limits. We estimate the CO excitation tempera- 
ture from the 2-1/1-0 ratio at selected outflow positions, which 
we find to range between 1 and 2.5. Since these ratios imply 
excitation temperatures between 7 and 20 K, we use a mean 
value of 15 K, as also suggested by the highest S/N CO(l-O) 
spectra, and which gives rise to the lowest column density for a 
given integ rated intensity. We als o assume a CO abundance of 
8.5 X 10-5 (lFrerkingetal.ll982h . 

The results of our outflow energetics are shown in Table 1 , 
and can be summarized by saying that the outflow has a mass 
of at least 0.01 Mq, a momentum of 0.1 Mq km s"', and an 
energy of 5 x 10^^ erg. A look at the distribution of these pa- 
rameters among the outflow components shows that most of 
the mass lies at low velocities while most of the momentum 
and energy is in the fastest gas. Quantitatively, the EHV com- 
ponent contains only 13% of the mass but carries 50% of the 
outflow momentum and more than 75% of its kinetic energy. 

Comparing the parameters of the blue and red gas (ratios 
in parenthesis in Table 1), we find a high degree of symme- 
try in the outflow energetics: the blue/red ratio is 1.2 for the 
EHV gas and even closer to 1.0 for the IHV component, de- 
spite the lower symmetry of its CO spatial distribution (Fig. 1). 
The SHV regime, on the other hand, has a factor of 3 more blue 
gas, but the study of this component is limited by the ambient 
emission. The energetics measurements, therefore, confirm the 
indications of the maps in Fig. 1 that at least the fastest part of 
the 104166 outflow has a very high degree of symmetry. 

To complete the description of the 104166 system we 
present in Fig. 3 NH^(l,l) and 1 .2mm c ontinuum maps of the 
IRAS source vicinity (see iMotte & Andre ■■2001. for a similar 
continuum map). These maps show that the dense core is nearly 
round and clearly peaked toward the IRAS source, and that it 
forms part of a SE-NW chain of cores (to be discussed in a 
future paper). Using standard NH3 analysis, we derive a gas 
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Fig. 3. Left: Integrated intensity map of the NH3 (1,1) emission 
toward the 104166 dense core. First contour and interval are 2 
K km s"' Right: Map of 1.2mm continuum emission toward 
the 104166 core (from Tafalla et al., in preparation) with sev- 
eral contours of the EHV outflow map from Fig. 1 . The contin- 
uum emission has been convolved with a 10" beam to enhance 
sensitivity, and its first contour and interval are 20 mJy per 11"- 
beam. Offset origin as in Fig. 1. 



kinetic temperature of 10 K, a turbulent component of 0.2 km 
s ' (FWHM), and assu ming a typical Taurus average para-NHs 
abundance of 5 x 10"^ JTafalla et al. I I2004S. we derive a mass 
of 1 .5 Mq, all typical of a Taurus core (Jiiina et al. 1999 ). 



4. Discussion 

The most striking element of the 104166 outflow is the EHV 
gas, with its high degree of collimation, large radial velocity, 
and strong blue/red symmetry. Because of its limited location 
and high velocity dispersion (> 10 km s"' , Fig. 2), the EHV gas 
is unlikely to represent the (still mysterious) agent that drives 
the outflow, although its extreme properties suggest that it must 
be closely connected to the outflow driving agent. Given the ge- 
ometry of the EHV gas (Fig. 1), this agent has to be jet-like (or 
have an important jet-like component), and in this context it 
seems natural to interpret the EHV gas as ambient cloud ma- 
terial prompt-entrained by a jet b ow shock. N umerical simu- 
lations of jet-driven outflows (e.g.. lLee et al. 112001.) . however, 
predict that gas accelerated by a bow shock will present a PV 
diagram with a "spur" feature having the fastest gas at the tip 
of the jet, in contrast with the almost constant velocity of the 
separate EHV component seen in Fig. 2. A possible solution 
to this discrepancy is the existence of multiple bow shocks, as 
suggested by the multiplicity of peaks in the maps of Fig. 1. 
High resolution observations of the 104166 outflow and realis- 
tic models of gas acceleration by time/angle variable jets are 
needed to test this interpretation and attempt to understand the 
origin of the EHV gas, a problem not peculi ar to the 10416 6 
system but common to all EHV outflows (e.g., Bachiller 1996). 
Comparing the 104166 outflow to other EHV outflows like 
L1448, 103282, or VLA1623 we conclude that it is a bone fide 
member of the class for its morphology, kinematics, and en- 
ergetics content. The 104166 outflow, in fact, seems a scaled- 
down version of the LI 448 outflow, as it contains about 5% of 



its momentum and energy while it has a central source with 
about 5% of its luminosity. This simple scaling (plus a similar 
kinematic age) shows that the 104166 outflow is almost as me- 
chanically efficient (L„,gc/Lt > 0.15) as L1448, and therefore 
belongs to the group of outflows an o rder of magnitude mor e 
powerful than typical Class 1 outflows ffiontemps et al. 192a)- 
The 104166 source, on the other hand, is a typi cal Taurus obj ect 
both in the parameters of its dust envelope iMotte & Andre I 
I2OOII) as well as in the turbulence and temperature of its core 
(Sect. 3.2). This implies that the (still unexplained) condi- 
tions required to produce an EHV outflow can be realized in 
the quiescent Taurus environment, and in addition, they can 
be achieved by a source with a luminosity as low as 0.4 Lq 
(the lowest of all EHV sources). Such a combination of cir- 
cumstances argues for EHV outflows resulting from a special 
(early) evolutionary stage in the protostellar phase, almost in- 
dependent of their luminosity and environment. 

Finally, we note that despite the similarities between the 
104166 outflow with the outflows from L1448, 103282, or 
VLA1623, the spectral energy distribution of 104166 is charac- 
terized by a bolometric temperature (Thoi) of 90-140 K, while 
the Thoi from the other sources with EHV outfl ows is closer to 
30-60 K JChen et al. I1 995HShirlev et al. l200nh . T his relatively 
large Ttoi yalue makes 104166 a Class I candidate JChen et al. I 
[1995), and therefore unexpectedly "evolved" compared with 
the other exciting sources of EHV flows. The recent non de- 
tection of a K band point source associated with 104166 by 
|Park & Kenvon (2002), however, suggests that 104166 may 
still be a truly Class source. A better characterization of 
104166 in the FIR with current (Spitzer) and future (Herschel) 
telescopes is therefore necessary. 
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